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ABSTRACT 

Context. Structural parameters of model star clusters are measured in radial profiles built from number-density, mass- 
density and surface-brightness distributions, assuming as well different photometric conditions. 

Aims. Determine how the core, half-star count and tidal radii, as well as the concentration parameter, all of which are 
derived from number-density profiles, relate to the equivalent radii measured in near-infrared surface-brightness and 
mass-density profiles. We also quantify changes in the resulting structural parameters due to depth-limited photometry. 
Methods. Star clusters of different ages, structure and mass functions are modelled by assuming that the radial distri- 
bution of stars follows a pre-defined analytical form. Near-infrared surface brightness and mass-density profiles result 
from mass-luminosity relations taken from a set of isochrones. Core, tidal and half-light, half-mass and half-star count 
radii, together with the concentration parameter, are measured in the three types of profiles, which are built under 
different photometric depths. 

Results. While surface-brightness profiles are almost insensitive to photometric depth, radii measured in number-density 
and mass-density profiles change significantly with it. Compared to radii derived with deep photometry, shallow pro- 
files result in lower values. This effect increases for younger ages. Radial profiles of clusters with a spatially-uniform 
mass function produce radii that do not depend on depth. With deep photometry, number-density profiles yield radii 
systematically larger than those derived from surface-brightness ones. 

Conclusions. In general, low-noise surface-brightness profiles result in uniform structural parameters that are essentially 
independent of photometric depth. For less-populous star clusters, those projected against dense fields and/or distant 
ones, which result in noisy surface-brightness profiles, this work provides a quantitative way to estimate the intrinsic 
radii by means of number-density profiles built with depth-limited photometry. 

Key words. Methods: miscellaneous; ( Galaxy.) globular clusters: general 



1. Introduction 

Star clusters are a powerful tool in the investigation of 
Galaxy structure and dynamics, star formation and evolu- 
tion processes, and as observational constraints to N-body 
codes. This applies especially to the long-lived and popu- 
lous globular clusters (GCs) that, because of their relatively 
compact nature, can be observed in most regions of the 
Galaxy, from near the center to the remote halo outskirts. 

In general terms, the structure of most star clusters 
can be described by a rather dense core and a sparse 
halo, but with a broad range in the concentration level. 
In this context, the standard picture of a GC assumes 
a isotherma l central region and a tid ally truncated outer 
region (e.g. iBinnev fc Merrifieldl fl998h . Old GCs, in par- 
ticular, can be v irtually considered as dy namically re- 
laxed systems (e.g. lNovola fe Gebhardtj[2~006D . During their 
lives clusters are continually affected by internal processes 
such as mass loss by stellar evolution, mass segregation 
and low-mass star evaporation, and external ones such as 
tidal stress and dynamical friction e.g. from the Galactic 
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bulge, disk and giant molecular clouds (e.g. Khalisi et alj 
2007; lLamers et alll2005t iGnedin fe Ostrikerlll997h . Over a 
Hubble time, these processes tend to decrease cluster mass, 
whic h may accelerate the core col lapse phase for some clus- 
ters (Di orgovski fe Mevlanl 11994 , and references therein) . 
Consequently, these processes, combined with the presence 
of a central black hole (in some cases) and physical condi- 
tions associated to the initial collapse, can affect the spatial 
distribution of ligh t (or mass) both in the central region and 
at lar ge radii (e.g. IGnedin et al.lll999i : iNovola fc Gebhardtl 
l2006h . 

It is clear from the above that crucial information re- 
lated to the early stages of Galaxy formation, and to 
the cluster dynamical evolution, may be imprinted in 
the present-day internal s tructure and large-scale spatial 
distribution of GCs (e.g. iMackev fc van den Berghl 12005; 
iBica et al.l[2006h . To some extent, this reasoning can be ex- 
tended to the open clusters (OCs), especially the young, 
which are important to determine the spiral arm and disk 
structures and the rotatio n curve of the Galaxy (e.g. iFriell 
119951 ; iBonatto et aill2006D . Consequently, the derivation of 
reliable structural parameters of star clusters, GCs in par- 
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ticular, is fundamental to better define their parameter 
space. This, in turn, may result in a deeper understand- 
ing of the formation and evolution processes of the star 
clusters themselves and the Galaxy. 

Three different approaches have been used to derive 
structural parameters of star clusters. The more traditional 
one is based on the surface-brightness profile (SBP), which 
considers the spatial distribution of the brightness of the 
component stars, usually measured in circular rings around 
the cluster center. The compilation of Harris (1996, and the 
2003 updat^H) presents a basically uniform set of parame- 
ters for 150 Galactic GCs. Among their structural param- 
eters, the core (R c ), half-light (RhL) and tidal (R t ) radii, 
as well as the concentration parameter c = log(Rt/R r ), 
were based mostly on the SBP database of iTrager et alj 
(|l995h . SBPs do not necessarily require cluster distances 
to be known, since the physically relevant information con- 
tained in them is essentially related to the relative bright- 
ness of the member stars. In principle, it is easy to mea- 
sure integrated light. However, SBPs are more efficient near 
the cluster center than in the outer parts, where noise and 
background starlight may be a major contributor. Another 
potential source of noise is the random presence of bright 
stars, either from the field or cluster members, especially 
outside the central region in the less-populous GCs or most 
of the OCs. Structural parameters derived from such SBPs 
would certainly be affected. One way to minimise this effect 
is the use of wide rings throughout the whole radius range, 
but this would cause spatial resolution degradation on the 
profiles, especially near the center. 

The obvious alternative to SBPs is to use star counts 
to build radial density profiles (RDPs), in which only the 
projected number-density of stars is taken into account, re- 
gardless of the individual star brightness. This technique 
is particularly appropriate for the outer parts, provided a 
statistically significant, and reasonably uniform, compari- 
son field is available to tackle the background contamina- 
tion. On the other hand, contrary to SBPs, RDPs are less 
efficient in central regions of populous clusters where the 
density of stars (crowding) may become exceedingly large. 
In such cases it may not be possible to resolve individual 
stars with the available technology. 

Finally, a more physically significant profile can be built 
by mapping the cluster's stellar mass distribution, which 
essentially determines the gravitational potential and drives 
most of the dynamical evolution. However, mass density 
profiles (MDPs) not only are affected by the same technical 
problems as the RDPs but, in addition, the cluster distance, 
age and a reliable mass-luminosity relation are necessary to 
build them. 

In principle, the three kinds of profiles are expected to 
yield different values for the structural parameters under 
similar photometric conditions, since each profile is sensi- 
tive to different cluster parameters, especially the age and 
dynamical state. Qualitatively, the following effects, basi- 
cally related to dynamical state, can be expected. Large- 
scale mass segregation drives preferentially low-mass stars 
towards large radii (while evaporation pushes part of these 
stars beyond the tidal radius, into the field), and high- mass 
stars towards the central parts of clusters. If the stellar 
mass distribution of an evolved cluster can be described 
by a spatially variable mass function (MF) flatter at the 
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cluster center than in the halo, the resulting RDP (and 
MDP) radii should be larger than SBP ones. The differ- 
ences should be more significant for the core than the half 
and tidal radii, since the core would contain, on average, 
stars more massive than the halo and especially near the 
tidal radius. Besides, the presence of brig ht stars preferen- 
tially in the central parts of young clusters (jBonatto fc Bical 
[20051 and references therein) should as well lead to smaller 
SBP core and half-light radii than the respective RDP ones. 

Another relevant issue is related to depth-limited pho- 
tometry. When applied to the observation of objects at dif- 
ferent distances, depth-limited photometry samples stars 
with different brightness (or mass), especially at the faint 
(or low-mass) end. Thus, it would be interesting to quan- 
tify the changes produced in the derived parameters when 
RDPs, MDPS and SBPs are built with depth-limited pho- 
tometry, as well as to check how the structural parameters 
derived from one type of profile relate to the equivalent 
radii measured in the other profiles. 

In the present work we face the above issues by deriving 
structural parameters of star clusters built under controlled 
conditions, in which the radial distribution of stars follows 
a pre-established analytical profile, and field stars are ab- 
sent. Effects introduced by mass segregation (simulated by 
a spatially variable mass function), age and structure are 
also considered. This work focuses on profiles built in the 
near-infrared range. The main goal of the present work is to 
examine relations among structural parameters measured 
in the different radial profiles, built under ideal conditions, 
especially noise-free photometry and as small as possible 
statistical uncertainties (using a large number of stars). In 
this sense, the results should be taken as upper-limits. 

This work is structured as follows. In Sect. [2] we present 
the star cluster models and build radial profiles with depth- 
limited photometry. In Sect. [3] we derive structural param- 
eters from each profile, discuss their dependence on depth, 
and compare similar radii derived from the different types 
of profiles. In Sect. [4] we compare relations derived from 
model parameters with those of the nearby GC NGC 6397. 
Concluding remarks are given in Sect. [5j 



2. The model star clusters 

For practical reasons, the model star clusters are simulated 
by first establishing the number-density radial distribution. 
The approach we follow is to build star clusters of differ- 
ent ages and concentration parameters, with the spatial 
distribution of stars truncated at the tidal radius (Rt). 
Stars are distributed with distances to the cluster center 
in the range < R < Rt, with the R coordinate h aving 
a num ber- frequency given by a function similar to a King 
(|1962f ) three-parameter surface-brightness profile. The mass 
and brightness of each star are subsequently computed ac- 
cording to a pre-defined mass function and mass-luminosity 
relation consistent with the model age. The last step is re- 
quired for the derivation of the MDP and SBPs. 

We point out that different, more sophisticated an- 
alytical models have also been used to fit the SBPs o f 
Galactic and extra-Galactic GCs, other than iKind fl962) 
profile. The most commonly used a re the single-mass, mod- 
ified isothermal sphere of iKind fl966) tha t is the ba- 
sis of the Galactic GC parameters given by ITrager et al. 
(|1995f ) and H03, the modified isothermal sphere of IWilson 
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Table 1. Model star cluster specifications 
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Table Notes. Col. 3: concentration parameter c = log(R t /R c ). Cols. 4 and 5: mass function slopes at the cluster center and tidal 
radius. Cols. 8-10: lower, upper and average star mass. Col. 11: absolute J magnitude at the turnoff (TO). Cols. 12 and 13: 
absolute J magnitude at the bright and faint ends. 



(|1975| ). that assumes a pre-defined stellar distribution func- 
tion (which results in more extended e nvelopes thanlKind 
119661) . and the power-law with a core o f lElson et all (11987D 
that has been fit to massive young clusters especially in 
the Magellanic Clouds (e.g. Mackey & Gilmore 2003a,b,c). 
Each function is characterised by different parameters that 
are somehow related to the cluster structure. However, the 
purpose here is not to establish a "best" fitting function of 
the structure of star clusters in general. Instead, we want to 
quantify changes in the structural parameters, derived from 
RDPs, MDPs and SBPs of star clusters with the stellar dis- 
tribution assumed to follow an analytical function, under 
different photometric conditions. We expect that changes 
in a given parameter should have a small dependence, if 
any at all, on the adopted functional form. 

The adopted King- like radial distribution function is ex- 
pressed as 



dN 



2nRdR 



y/1 + {R/R c ) 2 y/l + {Rt/Rcf 



(1) 



where ag is the projected number-density of stars at the 
cluster center, and R c and Rt are the core and tidal radii, 
respectively. Since structural differences are basically con- 
trolled by the ratio Rt/R c , we set R c = 1 in all models. Such 
a King- like RDP (for (T = 1-0) is obtained by numerically 
inverting the relation (see App. lA")) 
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where x = R/R c and u 2 = 1 + (Rt/R c ) 2 - Thus, a ran- 
dom selection of numbers in the range < n < 1 produces 
a King-like radial distribution of stars with the radial co- 
ordinate in the range < R/Rt < 1. The R/Rt curves as 
a function of n for the models considered in this work are 
shown in Fig. [1] (Panel a) . 

Once a given star has been assigned a radial coordinate, 
its mass is computed with a probability proportional to the 
mass function 



dN 
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where the slope varies with R according to x = x(R) = 
Xt + (Xt - Xo){R/Rt - 1), where xo and xt are the mass 



function slopes at the cluster center and tidal radius, re- 
spectively (Table [T]and Fig.[T]). Thus, the presence of large- 
scale mass segregation in a star cluster can be characterised 
by a slope xo flatter than xt- Mass values distributed ac- 
cording to Eq. [3] are obtained by randomly selecting num- 
bers in the range < n < 1 and using them in the relation 
of mass with n and x (App. [A| 

_ U, (m s /rrii) n , for x = 0.0, . . 

1 m s /[(l — n)(m s /m,i) x + n] 1 ^, otherwise, 

where m, and m s are the lower and upper mass values 
considered in the models (Table [T|). 

In what follows we adopt the 2MASS3 photometric sys- 
tem to build SBPs. Finally, the 2MASS J, H and K s mag- 
nitudes for each star are obtained according to the mass- 
luminosity relation taken fro m the correspondin g model 
(Table [J) Padova isochrone (|Girardi et al.l l2002ft . For il- 
lustrative purposes the model isochrones are displayed in 
Fig. [1] (panel c). 

The set of models considered here is intended to be ob- 
jectively representative of the star cluster parameter space. 
For globular clusters we use the standard age of lOGyr 
and the spatially uniform metallicity [Fe/H] = —1.5, which 
is typ ical of the metal-poor Galactic GCs (e.g. iBica et ail 
2006). However, we note that abundance variations have 
been suggested to occur within GCs (e.g. iGratton et alj 
2004). Basically, small to moderate metallicity gradients 
would produce slight changes in the colour and magni- 
tude of the stars in different parts of the cluster, which 
has no effect on the (star-count derived) RDPs and MDPs. 
The effect on the SBPs may be small as well, provided 
that the magnitude bin used to build the SBPs is wide 
enough to accommodate such magnitude changes. As for 
the core/tidal structure we consider the ratios Rt/R c = 
40, 20, 15, and 5, or equivalently the concentration pa- 
rameters c = log (Rt/Rc) ~ 1.6, 1.3, 1.2, and 0.7, which 
roughly correspond to the peaks in the distribution of c val- 
ues presented by the regular (non-post core collapse) GCs 
given in H03 (Fig. ED panel d). Models GC-A, B and D take 
into account mass segregation by means of a flat (xo = 0.00) 
mass function at the center and a Salpeter (1955) IMF 
(Xt — 1.35) at the tidal radius. GC-C model is similar to 
GC-B, except that it considers a uniform, heavily depleted 
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Fig. 1. Model star cluster specifications. Panel (a): a ran- 
dom selection of n in the range < n < 1 produces King- 
like RDPs in the range < R < R t (see Eq.|5J}. Panel (b): 
Radially-variable mass function slopes (4^- oc m~( 1+x ^ 
used in the models. Panel (c): Padova isochrones used to 
simulate the mass-luminosity relation of the star cluster 
models. The lOGyr, [Fe/H] = —1.5 metallicity isochrone is 
adopted in the globular cluster models. Panel (d): distribu- 
tion of concentration parameters of the GCs in H03 with 
peaks at c » 1.6, 1.3, and 0.7. Panel (e): model fraction of 
stars brighter than Mj — Mj TO + Ato- ln all cases, the 
fraction of stars brighter than the TO (Mj TO ) is below the 
1% level. 



MF (xo = 0.00) throughout the cluster. OCs are repre- 
sented by solar-metallicity models with the ages lOMyr 
(to allow for the presence of bright stars in young OCs), 
100 Myr (somewhat evolved OCs) and 1 Gyr (intermediate- 
age OCs), Rt/R c = 15 (c » 1.2) and a spatially variable 
MF (Table EJ. The values of c and the core/halo MF slopes 
are representative of OCs ()Bonatto fc Bica 2005) . Another 
effect not considered here is differential absorption. In prin- 
ciple, low to moderate differential absorption should have a 
minimum effect on the radial profiles, because of the same 
reasons as those given above for the metallicity gradient. 
High values, on the other hand, would affect RDPs as well, 
because of a radially-dependent loss of stars due to depth- 
limited photometry. However, inclusion of this effect is be- 
yond the scope of the present work. 

As expected, the fraction of stars brighter than the 
turnoff (TO) in the resulting star cluster models is signifi- 
cantly smaller than 1% (Fig. [TJ panel e). Thus, we had to 



use a total number of stars of 1 x 10 9 in all models, so that 
the radial profiles resulted statistically significant (small la 
Poisson error bars) especially at the shallowest magnitude 
depth. 

2.1. Depth-varying radial profiles 

The radial profiles were built considering all stars brighter 
than a given magnitude threshold, with the TO as refer- 
ence. At the bright end, statistically significant GC profiles 
were obtained for Ato = Mj,th—Mj,TO — —5, where Mjjh 
and Mj^to are the threshold and TO absolute magnitudes 
in the 2MASS J band. At the faint end, GC-models have 
Ato = 6.3. OC models have Ato = — 3 and — 4 at the 
bright end, and Ato = 8.9, 11.0, and 14.0, at the faint 
end. 

Starting at the bright magnitude end, RDPs, MDPs and 
SBPs were built considering stars with the J magnitude 
brighter than a given faint threshold, with the magnitude 
depth increasing in steps of Ato = 1 , up to the respective 
faint magnitude end. 

Figure O displays a selection of profiles corresponding 
to both extremes in magnitude depths, for the GC-D and 
OC-C models. These profiles are representative of the whole 
set of models, especially in terms of the small uncertainties 
associated with each radial coordinate. Reflecting the large 
differences in the number of stars at different photometric 
depths, the central values of the number and mass densities, 
and surface-brightness, vary significantly from the shallow- 
est to the deepest profiles. 



3. Structural parameters vs. photometry depth 

The depth-varying model SBPs are fit with t he em pirical 
three-parameter function introduced bv lKind (jl962l ) to de- 
scribe the surface-brightness distribution of GCs, which is 
characterised by the presence of the core and tidal radii. 
For RDPs and MDPs we use the King-like analytical pro- 
file that describes the projected number-density of stars as 
a function of R c and R t , cr(R) = ^FJTTr.' as gi ven by eq-ffl 



1-nRdRi 

We also compute the distances from the center which con- 
tains half of the cluster's total light, stars and mass. The 
half-star count (Rhsc), light (RhL) and mass (Rhm) radii 
are derived by directly integrating the corresponding pro- 
files. 

A selection of the resulting structural parameters as a 
function of Ato is given in Table O For simplicity we only 
present the values obtained from the bright and faint mag- 
nitude ranges, as well as for Mj < Mj.to- The whole set of 
parameters are contained in Figs.[3]-[H] At first glance, RDP 
and MDP radii present a significant decrease for shallower 
photometry, with respect to the intrinsic values. SBP radii, 
on the other hand, are more uniform. The most noticeable 
feature is that, except for GC-C (uniform mass function), 
RDP and MDP radii tend to become increasingly larger 
than SBP ones with increasing photometric depth. 



3.1. Dependence on photometric depth 

In Fig. [3] we compare the radii measured in GC profiles 
built with a given photometric depth (e.g. R c (Ato)) with 
the intrinsic ones, i.e. those derived from the deepest pro- 
files (R c ,deep)- RDP parameters are more affected than the 
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Table 2. Model star cluster structural parameters for different photometric depths 
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Model: GC-B; Input RDP parameters: R c = 1.0, Rt = 20.0 



-5.0 0.87 ±0.01 
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±6.3 1.00 ±(f) 
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Model: GC-C; Input RDP parameters: R c = 1.0, R t = 20.0 
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Model: GC-D; Input RDP parameters: R c = 1.0, R t = 40.0 



-5.0 0.90 ±0.01 
0.0 0.86 ±0.01 
±6.3 1.00 ±(f) 



2.81 ± 0.02 33.96 ±0.19 0.90 ± 0.01 2.81 ± 0.02 33.96 ±0.19 0.91 ± 0.01 2.82 ± (f) 34.18 ±0.05 

2.82 ±(f) 37.15 ±0.17 0.86 ± 0.01 2.82 ± (t) 37.15 ±0.17 0.91 ± 0.01 2.82 ± (f) 34.00 ± 0.05 
3.30 ±(f) 39.99 ±0.01 0.96 ± (t) 3.14 ± (t) 38.51 ±0.07 0.91 ± 0.01 2.82 ± (t) 34.20 ±0.04 



Model: OC-A; Input RDP parameters: R c = 1.0, R t = 15.0 

- : : — : — — = Q.81 ± 0.01 



-3.0 0.82 ±0.01 1.70 ±(f) 12.85 ±0.07 
0.0 0.78 ±0.01 1.72 ±(f) 13.78 ±0.06 
±8.9 1.00 ±(f) 2.08 ±(f) 



1.70 ±(f) 12.85 ±0.07 
1.72 ±(f) 13.78 ± 0.06 0.82 ± 0.01 
15.0 ±0.01 0.91 ±(f) 1.93 ±(f) 14.43 ±0.03 0.81 ± 0.01 
Model: OC-B; Input RDP parameters: R c = 1.0, R t = 15.0 

" 0.76 ±0.01 



0.82 ±0.01 
0.78 ±0.01 
0.91 ± (f) 



1.72 ±(f) 

1.72 ±(f) 

1.73 ±(f) 



13.18 ±0.02 

13.19 ±0.02 

13.20 ±0.01 



-3.0 0.72 ±0.01 
0.0 0.70 ± 0.02 
±11.0 1.00 ±(f) 



1.61 ± (f) 
1.61 ± (f) 
2.08 ± (f) 



13.30 ± 0.08 
13.67 ±0.08 
15.00 ±(f) 



0.72 ±0.01 
0.70 ±0.02 
0.84 ±0.01 



1.61 ±(f) 
1.61 ±(f) 
1.84 ±(f) 



13.30 ± 0.08 
13.67 ±0.08 
14.33 ± 0.04 



0.77 ± 0.01 
0.74 ± 0.02 



1.61 ±(f) 
1.61 ±(f) 
1.63 ±(f) 



12.75 ±0.03 

12.74 ±0.03 

12.75 ±0.02 



Model: OC-C; Input RDP parameters: R c = 1.0, R t = 15.0 



-4.0 0.62 ±0.02 
0.0 0.62 ± 0.02 
±14.0 1.00 ±(f) 



1.49 ±(f) 
1.49 ± (f) 
2.08 ± (f) 



13.06 ±0.10 
13.10 ±0.10 
15.00 ±(f) 



0.62 ±0.02 
0.62 ±0.02 
0.70 ± 0.02 



1.49 ±(f) 
1.49 ± (f) 
1.70 ±(f) 



13.05 ±0.10 
13.09 ±0.10 
14.35 ± 0.05 



0.71 ±0.01 
0.71 ±0.01 
0.64 ± 0.03 



1.49 ±(f) 
1.49 ± (f) 
1.49 ±(f) 



11.99 ±0.03 
11.99 ±0.03 
12.00 ±0.02 



Table Notes, (f): uncertainty smaller than 0.01 arbitrary units (au). The half-type radii are half-star counts (Rhsc)j half- mass 
(Rhivi) and half-light (RhL). 



MDP ones, while the SBP ones are essentially uniform, thus 
insensitive to photometric depth. Among the radii, RDP 
and MDP core are the most affected (underestimated) , fol- 
lowed by the half and tidal radii. In the most concentrated 
model (GC-A, c ss 0.7), measurements or R c in the RDP 
may be underestimated by a factor ss 25% in profiles shal- 
lower than near the TO, with respect to R c ,deep, and ~ 20% 
in MDPs. The effect is smaller in Rhsc an d RhM, which may 
be underestimated by « 15% in the same profiles. The un- 
derestimation in the tidal radii is smaller than ss 10%. As 
expected, RDP, MDP and SBP radii do not change when 
the mass function is uniform (GC-C model). 

Similar radii ratios in the OC models are examined in 
Fig. 3J Qualitatively, the same conclusions drawn from the 
GC models apply to the OC ones. However, the underesti- 
mation factor of RDP radii increases for younger ages, to 
the point that R c drops to w 60% of the deepest value for 
all profiles shallower than ss 3 mag below the TO in the 
OC-C model (lOMyr), and to w 70% for OC-B (100 Myr). 
The respective half-star count radii are affected by similar, 
although smaller, underestimation factors. MDP radii are 
less affected by cluster age than RDP ones. Similarly to 
the GC models (Fig. [3]), the three types of SBP radii are 



essentially insensitive to photometric depth, within uncer- 
tainties. We note that the presence of bright stars in the 
central region of young clusters (OC-C) appears to intro- 
duce a small dependence of the core radius on photometric 
depth (bottom-right panel). 

3.2. Comparison of similar radii among different profiles 

Differences on the same type of radii among the profiles, 
introduced essentially by a spatially variable MF, are dis- 
cussed in Fig.[5]for the GC models. Regardless of the model 
assumptions, RDP and MDP radii are essentially the same, 
except for the profiles corresponding to deep photometry, 
for which the RDP radii become slightly larger than the 
MDP ones. This occurs basically because of the larger frac- 
tion of low-mass stars at the outer parts of the clusters. 
Since all stars have equal weight in the building of the 
RDPs, the accumulation of low-mass stars at large radii 
ends up broadening the RDPs with respect to the MDPs. 
On the other hand, RDP core and half-star count radii 
tend to be larger than the SBP ones for profiles includ- 
ing stars fainter than near the TO. RDP R t may be 10 
- 20% larger than SBP ones for all depths. As discussed 
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Fig. 2. A selection of RDPs (top panels), MDPs (mid- 
dle) and 2MASS J magnitude SBPs (bottom) that illus- 
trate structural changes under different magnitude depths. 
Arbitrary units (au) are used both for the radial coordinate 
and projected area. 



above, the uniformly-depleted MF of GC-C model produces 
profiles whose radii are independent of photometry depth. 
The RDP to SBP core and tidal radii ratios decrease with 
concentration parameter. The RDP to SBP half-type radii 
ratios do not depend on c. 

The same analysis applied to the OC models is discussed 
in Fig. [6l The presence of massive stars in young clusters 
enhances the RDP to MDP radii ratios, especially the core 
and to some extent, the half- type radii. This occurs for pro- 
files that contain stars brighter than w 4 mag below the TO. 
For the youngest model (OC-C), the core radius measured 
in the RDP may be » 40% larger than the MDP one. This 
effect is enhanced when RDP radii are compared to SBP 
ones, again decreasing in intensity from the core to tidal 
radii. For OC-C, RDP core, half and tidal radii are w 55%, 
~ 40%, and « 25% larger than the equivalent SBP ones. 
Comparing with the GC models (Fig. [5]), the presence of 
a larger fraction of more massive (brighter) stars towards 
the center in young clusters tend to enhance radii ratios of 
RDP with respect to MDP, and especially, RDP to SBP. 

3.3. Further relations 

The models discussed in previous sections can be used as 
well to examine the dependence of the half-type radii with 
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Fig. 3. Structural parameters of the GC models. Top pan- 
els: Ratio of the tidal radius measured in profiles with a 
photometric depth Ato with respect to that derived from 
the deepest one, for the RDPs (left panels), MDPs (vertical- 
middle) and SBPs (right). Horizontal- middle panels: half- 
type radii. Bottom: core radii. TO values are indicated by 
the dotted line. Except for GC-C (uniform mass function), 
the remaining models present changes in radii in the RDPs 
and MDPs. SBP radii are essentially uniform. 



the concentration parameter, and to test how c varies with 
photometric depth. These issues are presented in Fig. [7] 

As already suggested by Figs. [3] and [J] the relation 
of the half radius with c, in a given model, changes sig- 
nificantly with photometric depth in RDPs (panel a) and 
MDPs (panel b). In SBPs, on the other hand, it is almost 
insensitive to depth (panel c). From eq. [TJ the half-star 
count radius is tightly related to the concentration param- 
eter according to R hS c = (0.69 ± 0.01) + (1.01 ± 0.01) c 2 . 
This curve fits well the values measured in the deepest RDP 
of all GC and OC models alike (panel a). Such a relation 
fails for the shallower profiles. A similar, but poorer, rela- 
tion applies to the values derived from the deepest MDPs 
(panel b), R hM = (0.63 ± 0.09) + (0.99 ± 0.05) c 2 . It fails 
especially for the young (OC) models. The GC SBPs, on 
the other hand, can be poorly fit with the linear function 
RhL = (-0.9 ± 0.1) + (2.4 ± 0.1) c (panel c). 

Concentration parameters measured in RDPs and 
MDPs (panels d and e) change with photometric depth. 
Around the TO they reach the maximum value, which cor- 
responds to a star cluster ss 15% more concentrated than 
the pre-established value (Table [1]). At the shallowest pro- 
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-4 4 8 12 -4 4 8 12 -4 4 8 12 

A T0 (mag) A TO (mag) A T0 (mag) 

Fig. 4. Same as Fig. [3] for the OC models. For comparison 
purposes, the y-scale is the same as in Fig. [3] Similarly to 
the GC models, radii changes are conspicuous in the RDPs 
and MDPs. 



files c presents a value intermediate between the maximum 
and the pre-established one, which is retrieved at the deep- 
est profiles with the inclusion of the numerous low-mass 
stars. The exception again is the uniform MF model GC-C, 
whose c values do not change with Ato- c values measured 
in SBPs are essentially insensitive to photometric depth 
(panel f). 

4. NGC6397: a test case 

We compare the results derived for the model star clus- 
ters with similar parameters measured in the My = —6.63, 
nearby GC (d Q = 2.3 kpc) NGC6397. Being populous is 
important to produce statistically significant radial profiles, 
while the proximity allows a few magnitudes fainter than 
the giant branch to be reached with depth-limited photom- 
etry. 

NGC 6397 is a post-core collapse GC with evidence of 
large-scale mass segregation, as indicated by a mass func- 
tion flatter at the center than outwards (|Andreuzzi et all 
2004 and references therein) . 

Additional relevant data (from H03) for the metal-poor 
([Fe/H] = -1.95) GC NGC 6397 are the Galactocentric dis- 
tance Rqc — 6 kpc, half-light and tidal radii (measured in 
the V band) R hL = 2.33' and R t = 15.81', and Galactic 
coordinates I = 338.17°, b = -11.96°. Thus, bulge star 
contamination is not heavy, and cluster sequences can be 
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Fig. 5. GC model profiles. Ratio between the same type 
of radii as measured in RDPs and MDPs (left panels) and 
RDPs and SBPs (right panels). From top to bottom: tidal, 
half and core radii. TO values are indicated by the dotted 
line. 

unambiguously detected, which is important for the extrac- 
tion of radial profiles with small errors (see b elow). Using 
SBPs b uilt with 2MASS im ages and a fit with|Kingl (|1962|) 
profile, ICohen et al.l (|2006h derived the core radius in the 
J band R C (J) = 61.5" ± 9.3". However, based on Hubble 
Space Te lescope data and using a po wer-law plus core as fit 
function. iNovola & Gebhardtl (|2006f ) derived R c = 3.7" in 
the equivalent V band, thus roughly resolving the post-core 
collapse nucleus. 

The post-core collapse state of NGC 6397 does not affect 
the present analysis, since the goal here is the determina- 
tion of changes produced in cluster radii derived under the 
assumption of a King-like profile (Sect. [3]) applied to RDP, 
MDP and SBPs built with different magnitude depths. We 
base the analysis of NGC 6397 on J, H and K s 2MASS pho- 
tometry extracted using VizieR[f| in a circular field of radius 
Rcxt = 70' centered on the coordinates provided in H03. 
This extraction radius is large enough to encompass the 
whole cluster, allowing as well for a significant comparison 
field. 

For a better definition of the cluster sequences we ap- 
ply the statistic a l deco ntamination algorithm described in 
iBonatto fc Bica! (|2007f ). which takes into account the rela- 
tive number-densities of candidate cluster and field stars in 
small cubic CMD cells with axes corresponding to the mag- 



vizier.u-strasbg.fr/viz-bin/VizieR?-source=II/246 
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Fig. 6. Same as Fig. [5] for the OC models. For comparison 
purposes, the y-scale is the same as in Fig. [5j 



nitude J and the colours (J — H) and (J — K s ). Basically, 
the algorithm (i) divides the full range of magnitude and 
colours of the CMD into a 3D grid, (ii) computes the ex- 
pected number-density of field stars in each cell based on 
the number of comparison field stars with magnitude and 
colours compatible with those in the cell, and (Hi) sub- 
tracts the expected number of field stars from each cell. 
Typical cell dimensions are AJ = 0.5, and A(J — H) = 
A(J — K s ) = 0.25, which are large enough to allow sufficient 
star-count statistics in individual cells and small enough 
to preserve the morphology of the CMD evolutionary se- 
quences. The comparison field is the region located between 
50 < R(') < 70, which is beyond the tidal radius. 

Field-decontaminated CMDs allow for a better defi- 
nition of colour-magnitude filters, useful to remove stars 
(and artifacts) with colours compatible with those of the 
field which, in turn, improves the cluster/background con- 
trast in RDPs and SBPs. They are wide enough to ac- 
commodate cluster MS and evolved star colour distribu- 
tions and dynamical evolution-related effects, such as en- 
hanced fractions of binaries and other multi ple systems (e.g. 
iBonatto fc Bicall2007t iBonatto et al.ll2005l ). 

Figured] (panel a) displays the decontaminated CMD of 
a central region of NGC 6397, with R < 5', somewhat larger 
than the half-light radius (Table[3]). We take J = 15 as refer- 
ence to extract the depth- variable profiles. RDPs and SBPs 
are built with colour-magnitude filtered photometry, with 
the faint end varying in steps of Ajis = 0.5, with the deep- 
est (i.e. at the available 2MASS depth) profile beginning at 



-5 5 10 -5 5 10 -5 5 10 
A TO (mag) A TO (mag) A TO (mag) 

Fig. 7. Top panels: relation of the half-type radii with the 
concentration parameter, for the RDPs (left panel), MDPs 
(middle) and SBPs (right). For each model, Rh values in- 
crease for deeper profiles. Dashed line in panels (a) and (b): 
Rh ~ c 2 . In panel (c): RhL ~ c. Bottom panels: concentra- 
tion parameter as a function of photometric depth. 



J = 16 and the brightest one ending near the giant clump at 
J = 13. The extracted profiles are fitted with the King-like 
function discussed in (Sect. [2]). A selection of depth- limited 
RDPs, together with the respective fits and uncertainties, 
is shown in Fig. and the corresponding RDP and SBP 
(J band) radii are given in Table [3J Within uncertainties, 
the present value of the core radius (for the dee pest profile), 
R C (J) = 1.4' ±0.3', agrees with that derived by I Cohen et alj 
(|2006f ). using the same fit function. The near- infrared half- 
light radius, on the other hand, is larger than the optical 
one (H03), R hL (J) w 1.5R hL (^). 

Effects of the varying magnitude depth on the radii of 
NGC 6397 are examined in Fig. [5] Qualitatively, the result- 
ing curves agree, within uncertainties, with the behaviour 
predicted by the GC models (Figs. [3] and [5]). Compared to 
the values measured in the deepest RDP, the tidal (panel 
a), half-star counts (b) and core (c) radii decrease for shal- 
lower profiles, especially for Ajis > —0.5, remaining almost 
uniform for Ajis < —0.5. In particular, the core radius 
measured in shallow RDPs (containing essentially giants) 
drops to 45% of its deepest value (which includes stars 
at the top of the MS). Consistently with the GC models 
containing a spatially variable MF (Sect. [3]), the varying 
depth affects the tidal, half and core radii, with increasing 
intensity. SBP radii, on the other hand, remain essentially 
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Fig. 8. Structural analysis of NGC6397. Panel (a): de- 
contaminated CMD of a central (R < 5') region. The 
reference magnitude J = 15 is indicated by the dashcd- 
line. Shaded region: colour-magnitude filter. Background- 
subtracted RDPs for stars brighter than J < 15 + Aji5, 
with A,/i5 = 1, 0, — 1, — 2 are shown in panels (b) to 
(d), respectively. The respective King-like fits (solid- line) 
together with the fit uncertainty (shaded region) are shown. 



uniform with variable depth, consistent with the GC mod- 
els (Sect. [3]). The same conclusions apply to the RDP to 
SBP radii ratio (right panels). 



5. Concluding remarks 

In this work we simulated star clusters of different ages, 
structure and mass functions, assuming that the spatial dis- 
tribu t ion of stars follows an analytical function, similar to 
iKind (|1962f) profile. The mass and near-infrared luminosi- 
ties of each star were assigned according to a mass func- 
tion with a slope that may depend on distance to cluster 
center. They form the set of models from which we built 
number-density, mass-density and surface-brightness pro- 
files, allowing for a variable photometric depth. The struc- 
tural parameters core, half-light, half-mass and half-star 
count, and tidal radii, together with the concentration pa- 
rameter, were measured in the resulting radial profiles. Next 
we examined relations among similar parameters measured 
in different profiles, and determined how each parameter de- 
pends on photometric depth. We point out that the results 
should be taken as upper-limits, especially for open clus- 
ters, since we have considered noise-free photometry and 
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Fig. 9. Left panels: RDP and SBP structural radii of 
NGC6397 as a function of Ajis, normalised to the val- 
ues measured in the deepest profile. Right panels: RDP to 
SBP radii ratios (similar to Fig.[5|). 

Table 3. Radii of NGC 6397 from RDPs and 2MASS SBPs 
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Table Notes. 



( 1962) functions (Sect. [3]) to the respective profiles. The half- 
star counts and half-light radii were measured directly on the 
profiles. 



a large number of stars, which produced small statistical 
uncertainties. 

With respect to the adopted form of the radial dis- 
tribution of stars, we note that IKind (jl962t ) isothermal 
sphere, single-mass profile h as been supe r seded b y more 
real istic models like th ose of IKind (|1966f ). IWilsonl {l975) 
andU son et all (jl987f ). which have been fit mostly to the 
SBPs of Galactic and extra-Galactic GCs (Sect. [2}. The 
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analytical functions associated with these models are char- 
acterised by different scale r adii (among other parameters) 
that are roughly related to iKind |l962) radii. Thus, it is 
natural to extend the s caling with photometric depth un- 
dergone by [King (1962) radii to the equivalent ones in the 
other models. 

The main results can be summarised as follows. 

— (i) Structural parameters derived from surface- 
brightness profiles are essentially insensitive to photo- 
metric depth, except perhaps the cluster radius in very 
young clusters. 

— (it) Uniform mass functions also result in structural pa- 
rameters insensitive to photometric depth. 

— (Hi) Number-density and mass-density profiles built 
with shallow photometry result in underestimated radii, 
with respect to the values obtained with deep photome- 
try. Tidal, half-star count and half-mass, together with 
the core radii are affected with increasing intensity. 

— (iv) Because of the presence of bright stars, radii under- 
estimation increases for young ages. 

— (v) For clusters older than ~ 1 Gyr, number-density and 
mass-density radii present essentially the same values; 
for younger ages, RDP radii become increasingly larger 
than MDP ones, especially at the deepest profiles. 

— (vi) Irrespective of age, profiles deeper than the turnoff 
have RDP radii systematically larger than SBP ones, 
especially the core. 

— (vii) The concentration parameter also changes with 
photometric depth, reaching a maximum around the 
turnoff. 

Most of the above model predictions were qualitatively 
confirmed with radii measured in ground-based RDPs and 
SBPs of the nearby GC NGC 6397. 

In principle, working with SBPs has the advantage of 
producing more uniform structural parameters, since they 
are almost insensitive to photometric depth. However, as 
discussed in Sect. [U SBPs usually present high levels of 
noise at large radii. Noise that is also present in SBPs of 
clusters projected against dense fields and/or the less pop- 
ulous ones. A natural extension of this work would be to 
examine radial profiles built with photometry that includes 
observational uncertainties, differential absorption, metal- 
licity gradients, binaries, and star cluster models with a 
number of stars compatible with those of open clusters. 

As a consequence of the wide range of distances to the 
Galactic (and especially extra-Galactic) star clusters, in- 
terstellar absorption, and intrinsic instrumental limitations, 
the available photometric data for most clusters do not sam- 
ple the low-mass stars. All sky surveys like 2MASS, usually 
are restricted to the giant branch, or the upper main se- 
quence, for clusters more distant than a few kpc. In such 
cases, the structural parameters have to be derived from 
radial profiles built with photometry that does not reach 
low-mass stars. The present work provides a quantitative 
way to estimate the intrinsic (i.e. in the case of photome- 
try including the lower main sequence) values of structural 
radii of star clusters observed with depth-limited photom- 
etry. 
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Appendix A: Transformation of a random number 
distribution into an analytical function. 

The simulations discussed in the present work depend on 
the transformation of a distribution of numbers p(n) , with n 
varying in the range < n < 1 , into the analytical function 
f{x), with x in the range x min < x < x max . What results 
from this is that a random sorting of numbers n would 
produce x variables with a number-frequency given by the 
pre-defined function f(x). Formally, the transformation is 
expressed as 

f(x)dx—p(n)dn. (A.l) 

In a random distribution, all numbers have the same 
probability, so that we can set p(n) — 1. Thus, the formal 
relation of n to x is given by 

n(x)= f(x')dx'/ f(x')dx'. (A.2) 

J Smin * ^min 

The analytical (or numerical) inversion of the latter re- 
lation gives x = x(n). 

For example, consider a mass function <j)(m) — cx 
mr( 1+x \ with the mass varying in the range rrii < 
w(M Q ) < m s . If x 7^ 0.0, what results from Eg. IA.2I is 

n(m) — (m~' x — m~ x ) / (m~ x — m^ x ) . 
Solving this for m produces 

m(n) = m s /[(l — n)(m s /mi) x + n] l ^ x . 
For x — 0.0 the result is m(n) = rrii (m s /mi) n . 
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In the case of the three-parameter King-like profile 
(Eq. [T]in Sect. [5]), the first step in the transformation, rep- 
resented by Eq. IA.21 results in the n = n(R) function given 
by Eq. [2] (Sect. [2]). However, the latter relation cannot be 
analytically inverted, so that this task has to be done nu- 
merically, using the curves plotted in Fig. Q] (panel a). 



